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O ■ Abstract: We study the bound states of brane/antibrane systems by examining the 

motion of a probe antibrane moving in the background fields of source branes. The 
p • classical system resembles the point-particle central force problem, and the orbits can 

D . be solved by quadrature. Generically the antibrane has orbits which are not closed on 

themselves. An important special case occurs for some Dp-branes moving in three 



transverse dimensions, in which case the orbits may be obtained in closed form, 
giving the standard conic sections but with a nonstandard time evolution along the 
orbit. Somewhat surprisingly, in this case the resulting elliptical orbits are exact 
solutions, and do not simply apply in the limit of asymptotically-large separation 
or non-relativistic velocities. The orbits eventually decay through the radiation of 
massless modes into the bulk and onto the branes, and we estimate this decay time. 
Applications of these orbits to cosmology are discussed in a companion paper. 
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1. Introduction 

One of the most important developments in string theory of late has been the re- 
alization that it contains many kinds of extended objects like Dp-branes. In many 
situations these objects can be treated much like elementary particles, in that they 



- 1 - 



carry a generalised mass and charge and these give rise to mutual forces which govern 
their relative motion. 

This observation is all the more interesting when the dimension of the brane 
exceeds three, since then the moving branes can contain an entire world inside of 
themselves (possibly including us). In this kind of picture the collision of Dp-branes 
can be interpreted as colliding universes. This observation has recently been used 
to propose several novel scenarios for the very-early universe, giving rise to inflation 
or alternatives to inflation ||, ||, |^, 0, |, ^, allowing in both cases an at- 
tractive geometrical and stringy interpretation of the various cosmologically-active 
scalar fields. 

Perhaps the best-studied brane interactions are brane collisions, which in the 
particular case of branes and antibranes can lead to brane annihilation (for which 
insights have come from recent studies of tachyon condensation). However these 
kinds of direct collisions only begin to explore the rich kinds of physics which are 
possible when particles and branes mutually interact. A different — and for Dp- 
branes as yet largely unexplored — possibility is the development of bound orbits, 
such as commonly arise for particles interacting through central forces. 

In this paper we begin the study of the properties of the brane generalisation 
of such bound orbits. We do so partly motivated by the potential applications of 
such a study to brane-world cosmology. In particular we could imagine the brane 
which makes up our own universe being in a closed orbit around another brane 
due to the mutually attractive forces which act between them. But regardless of 
these applications, we regard brane bound states to be of interest in their own right 
inasmuch as they add to our understanding of how branes and strings can interact. 

We start this study here with a discussion of the bound states which are possible 
for brane/antibrane systems. These are among the simplest to use in applications, 
partly because of the triviality of the various consistency conditions which would 
otherwise arise when trying to construct brane configurations inside compact extra 
dimensions. We call this kind of bound state 'branonium' due to its obvious similarity 
with the analogous electron/positron system: positronium. 

Our discussion is couched in terms of a test antibrane which is moving within 
the fields sourced by a system of parallel branes. We set up the effective brane 
lagrangian for this motion and analyse the features of the resulting bound orbits. It 
happens that despite the many complications of the resulting effective Lagrangian, 
it is still possible to solve explicitly for the trajectory of the extended test brane, 
with a result which is very similar to those obtained for nonrelativistic motion in a 
central force. In particular we find the conditions for the existence of stable closed 
orbits, and show that these include nontrivial systems of D-branes, with perfectly 
elliptical (non-precessing) orbits being possible if the branes move in three transverse 
dimensions. 

Finally, we give a preliminary study of the stability of the system due to different 
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sources, including radiation reaction due to radiation of bulk and brane modes, as 
well as nonperturbative instabilities and tachyonic modes. The dominant instability 
appears to be towards radiation into massless modes, and we compute the orbital 
decay time for this process. 

2. Probe-Brane Analysis 

We start by setting up and solving the equations which are relevant to the motion 
of a probe antibrane moving within the background of N parallel branes. To do so 
we first examine the fields which the source branes set up, and then examine the 
equations of motion which they imply for a test brane or antibrane. We show how 
these equations may be solved for the motion of the test-brane centre of mass when 
all other brane modes are frozen. (We return to the question of why these modes 
should be frozen in the next section.) Finally we examine the bound-state orbits and 
identify the surprising special class for which the exact relativistic orbits are ordinary 
conic sections, and we study the properties of the orbits in this special case in some 
detail. 

2.1 The Fields Due to Source Branes 

Let us first consider the gravitational, dilaton and 'electromagnetic' fields which are 
set up by a set of N parallel p-branes in L> -dimensional spacetime. We take the 
action for these fields to be given in the Einstein Frame by 

I ^""^ {Rmn + 9m 9^ 0) + 2^ e"^ Fm,...m„ F^-'^r. ^ 

(2.1) 

where the n-form field strength is related to its {n — l)-form gauge potential in the 
usual way F^n] = dA[n-i]- Here n is related to the spacetime dimension d of the N 
p-branes by d = p+ l = n — 1. We denote the d coordinates parallel to the branes 
by and the D — d transverse coordinates by y™- . 

The constant a depends on which kind of brane is being considered. For instance, 
if were to arise from a string-frame action 

with the string-frame metric given by Qmn = ^^"^ Qmn for A = 4/(D — 2), then 
a = as + 2 — \{n — 1) . Two cases of particular interest are: {i) NS-NS fields (like the 
fields rising in the gravity supcrmultiplet in various dimensions), for which Us = —2 
and so a = aNS = — 4(n — 1)/{D — 2); {ii) R-R fields, for which = and so 
a — aR — 2{D — 2n)/{D — 2). In the special case D = 10 we have A = | and so 
these two cases become aNs — — |(^~1) — ~ |(p + l) or oiR — 1(5 — ?^) = |(3— p), 
respectively. 
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The solutions to the equations for this system which describe the fields due to a 
stack of parallel collocated p-branes is (in the Einstein Frame) |T0[| : 



ds^ = h-Hx^ + h^dy^ , = h\ Am,„p = C(l - h'^) , (2.3) 

where all the other components of the {n — l)-form field vanish. The constants 7, 
7, K and ^ are given by 

4d 4d 2a ^ 2 , ^ 

7 = 77^^TVX' 7= ,n ON A ' '^ = ^' C = ^, (2.4) 



(D-2)A' ' (D-2)A' A' ^ 

where d = D — d — 2 and 

. r> 2 , . 

We assume throughout the existence of more than 2 transverse dimensions, so d > . 
The quantity h denotes the harmonic function: 

h{r) = l + k/r^ , (2.6) 

where r is the radial coordinate transverse to the branes, = 5mn y"^y"' , and k is 
an integration constant. 

Dp-branes in 10 dimensions represent a particularly interesting special case of 
this solution, corresponding |Tl| to the choices a = = ^{3~p) , A = 4, d = p + 1, 



d = 7 — p, K = i(3 — p) and k = c^QsN Ig'^ . (Clearly the condition c? > then 
implies p < 7.) Here Is is the string length, which is related to the string tension 
by Is^ = a' . The string coupling constant at infinity is denoted gs and Cp is the 
constant Cp = (20?)^-^ r(^) . 

Having an explicit string solution is useful for many reasons, not the least of 
which is the understanding which this gives of the limits of validity of the field 
theoretic description we shall use throughout what follows. For these purposes there 
are two dimensionless quantities which are relevant. The first of these controls the 
validity of string perturbation theory, which for string self-interactions requires the 
local string coupling, g{r) = QsG'^, to be everywhere small. The smallness of this 
quantity is required in order to suppress the modifications which higher string loops 
would otherwise have on our semiclassical analysis. 

The second important dimensionless quantity is the ratio r/lg, which must be 
large in order to believe that the brane interactions are dominated only by the me- 
diation of bulk fields which correspond to massless string states. Once r ~ /<, higher 
order corrections in a' become important because over these short distances contri- 
butions from higher string modes become possible. A particularly dangerous example 
of such a correction arises if a probe antibrane is introduced, since then there is an 
open string tachyon which appears at short distances within the open-string sector 
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for strings which stretch between the branes and the antibrane. Within the closed 
string picture this tachyon can be regarded as a kind of a' correction, since by dual- 
ity it can be thought of as emerging when all the closed string modes are resummed. 
The appearance of this tachyon signals the onset of the instability in which the probe 
antibrane annihilates with one of the source branes. 

One might worry that the requirement r ^ Ig might also preclude using the full 
nonlinear solution given above, requiring us instead to use only its leading approxima- 
tion in powers of k/r'^ . But this need not be required if the solution's characteristic 
length scale, p = k^/'^ = ls{cp gsNY/^'^~P\ is also large, p Z^, since in this case 
we may have k/r'^ = 0(1) even if r ^ Ig. Notice that for Dp-branes the condition 
p ^ Is requires gsN ^ 1 , and so takes us beyond the limit of open-string pertur- 
bation theory, which requires ggN -C 1 (where the factor of here arises from the 
Chan-Paton contributions from each world sheet boundary). In this way the field- 
theory limit provides (for r ^ Ig) results which are nonperturbative in the coupling 
of the bulk fields to the branes. 

2.2 The Probe Brane 

In the above background let us now follow the motion of another single parallel 
p-brane, displaced from the original by the radial coordinate-distance r (for a 
related discussion see ||14[)- This motion is described by the brane action, which 
can be decomposed as the sum of two pieces: the Born-Infeld (BI) and the Wess- 
Zumino (WZ) parts. The dilaton and graviton couplings are given by the Born-Infeld 
contribution, which is (in the String Frame): 

Sbi = -Tp J rft e-<^i/- det(^Af^ d^x^'d.x^ + 2 h^J^^,) , (2.7) 

where are the coordinates of the embedding, ^'^ are the world- volume coordinates 
and Tp denotes the brane tension. J^f^i, denotes the field strength for any open-string 
gauge modes confined to the brane, although for the moment we put these fields to 
zero. 

The coupling to the bulk gauge field, A^+i] , is given by the Wess-Zumino part 
of the brane action: 

Swz = -qTp J A[p^i], (2.8) 

where q represents the brane charge, with q = 1 representing a probe brane and 
q = —1 representing a probe antibrane. 

Since our branes are straight and parallel it is convenient to choose the following 
coordinate gauge: 

x'' = e, (2.9) 

where as before are the coordinates parallel to the world volume. To the extent 
that we follow only the brane's overall centre-of-mass motion — we return to the 
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validity of this assumption in the next section — we may also take 



(2.10) 



where t = denotes coordinate time and y"^{t) are the coordinates perpendicular to 
the brane surface, which parameterise the position of the probe. With these choices 
the induced String-Frame metric on the probe becomes: 

ds^ = e^^ ds^ = h-^ [{-1 + h^v^) de + [dCf] (2.11) 
where i = 1, and = '^rni^V"^ / ■ '^^^ constants /3 and u are given by 

4:{d-2a) 



/3 = 7 — kA 



UJ = J + J 



4 
A 



(2.12) 



(D-2) A ' 

For the special case of Dp-branes in 10 dimensions these become P = \ and u = 1. 



The effective action for the probe brane is then: 

1 



5* — TpVp 



dt 



h 



(2.13) 



where Vp = J d^^ is the spatial volume spanned by the probe brane, ( 
by eq. ( p. 41) and 

f3d _2 [a{d- d) +dd] 
''"''+1"" (D-2)A ■ 
For the special case of Dp-branes in 10 dimensions we therefore have rj = 1. 
Apart from an additive constant, the brane action has the form 

1 



uj given 



(2.14) 



S = —m / dt 



Vl - h^v"^ 



(2.15) 



with q = C q and m = TpVp . This is very similar to that of a particle moving in the 
D — d dimensions transverse to the branes. That part of the action containing the 
square root corresponds to the particle's relativistic kinetic energy and the coupling 
to the dilaton field, while the part involving q corresponds to the potential due to 
the Ramond-Ramond couplings. Notice that the apparent mass, m, of the particle 
scales with the brane world-volume, and so becomes infinitely large in the limit that 
this volume diverges. 

Several features of the probe-brane action bear emphasis: 

• A measure of the potential energy of the probe brane may be obtained by 
evaluating the brane action for a brane at rest: v = 0. The result obtained 
in this way is proportional to the combination /i^'' — q/h. For Dp-branes in 
10 dimensions this becomes {1 — q)/h, which is nonzero for a probe antibrane 
(g = —1), but vanishes for probe brane (g = 1). This is as one expects given 
the BPS nature of the Dp-brane. 
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• The brane motion is along a null geodesic in the transverse dimensions when 

= /i""^. This becomes v = 1 asymptotically far from the set of N source 
branes, but is smaller than unity closer to the source branes. 

• Masslcss particles trapped on the probe branc move along null geodesies of the 
induced metric, which corresponds to (d^'^ / dt)'^ = 1 — h^v"^ , indicating that the 
effective maximum speed for trapped particles, = 1 — h^v"^ , varies with time 
as the branc moves through the fields of the source branes. Notice that — >■ 
if the brane becomes ultra- relativistic: f ^ — > h~'^ . 

2.3 Classical Dynamics 

Let us analyse the movement va. — D — d dimensions of a particle with the 
Lagrangian: 

L = -m (^^Vl -h'^v'^- 1^ , (2.16) 



with h = 1 + k/r'^ and A; > 0. 

Invariance of this Lagrangian under 0{dT) rotations ensures the conservation 
of the angular momentum tensor, Lij = XiPj — XjPi, with Xi and pj denoting 
the components of the particle's position and canonical momentum vectors. The 
conservation of this tensor ensures that the particle motion is confined to the plane 
which is spanned by the particle's initial position and momentum vectors. Let us 
denote by r and 9 polar coordinates in the plane of the particle motion. The 
Lagrangian for the resulting two-dimensional system is 



L = -m (h-i^l - h"" (r2 + r^O^) - qh-^j . (2.17) 
Invariance under translations of 9 ensures the conservation of angular momentum 
^ = ^ = , (2.18) 

Notice that this is just the relativistic angular momentum. In terms of this the 
angular velocity becomes 



r 



The canonical momentum associated with the radial motion is similarly 

Pr r 



Pr - — 

m 



(2.20) 
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Inverting this to obtain r as a function of pr then gives 

Constructing the canonical Hamiltonian gives the energy formula: 



E _ 1 
m 



1 + h^"^' 

which is conserved during the motion. 



1/2 



1 
h 



(2.21) 



(2.22) 



2.3.1 Qualitative Analysis 

A qualitative understanding of the orbits which are implied by this system may be 
obtained using the effective potential 

1/2 



Vcs{r) = e{pr 



1 + /i'"-" -T 



1 
h 



(2.23) 



which we plot in Figure (|T|) for the case p = 6. 

The utility of this potential follows from the observation that £ is a monotonically 
increasing function of p,-'- de/dpr > 0. The allowed range of r for classical motion 
is therefore easily found by the standard device of plotting Vcs{r) against r, and 
finding those r's for which e > Vcs{r) . 

The behaviour of Vcs follows from the properties of h, which is a monotonically 



decreasing function of r, with limits h 
r ^ oo. This implies 



k/r" 



CX3 as r 



and h 



eflf 



(1-^) + 



(g — ifjk/r'^ 



for c/ < 2 
for J = 2 
for > 2 



as r 



oo 



1 as 



(2.24) 



V, 



eff 



for {uj-2)d<2 
for (cu - 2) d = 2 
for {uj-2)d>2 



as r 



0. 



(2.25) 



and, provided rj and uj are positive, 

[E/k^/^-q/ky 
—qr'^jk 

These limiting forms are often sufficient to determine the existence of a minimum 
for \4fj for some value r = r^, which signals the existence of a stable circular orbit 
at this radius. For instance, for the 10 dimensional Dp-brane example we have 
ijj = ri = \ and (i = 7 — p > . This then gives the limits 

{((/ — l)k/r for p = 6 

[(."^ /2 + q-l\k/r'^ for p = 5 as r ^ oo , 

£V(2r2) for p < 5 
Kff ^ Ir^^-'p'^/^/k^''^ as r ^ 0. (2.26) 
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These expressions show that 
for p > 5 Ves blows up (to pos- 
itive infinity) at the origin, and 
goes asymptotically to a constant 
as r ^ oo. Furthermore, if g = — 1 
then Ves begins decreasing as r 
comes in from infinity, and so there 
must be a minimum somewhere 
in between. The position, r^, of 
the minimum is the solution to the 
equation (2rc + 1)^^^ = 8r^. 

For p < 5, Ves instead ap- 
proaches zero both as r ^ and as 
r ^ oo , and it grows as r comes in 
from infinity. The effective poten- 
tial in this case instead has a max- 
imum at some intermediate radius, 
and no stable circular orbits exist. 
Bound motion exists for distances 
smaller than the one corresponding 
to the maximum of the potentials 
but the trajectories are not closed 
and eventually pass through the 
origin. 



4- 



3- 



V(r) 



Effective potential for D6 and antiPS-brane probe 




10 



Legend 



antiD6-brane probe 
D6-brane probe 



Figure 1: The effective potential for the radial 
motion of the brane centre-of-mass for two differ- 
ent probes, a D6 and an anti-D6-brane in the back- 
ground of a large number N of D6-branes. We 
have taken Ngjg = 2. Bound states exist only for 
anti-D6-branes. 



The angular frequency, Q, for a circular orbit is a constant, which is related to 



the orbital radius, 



by eq. ( p.l9| ) specialized to r 



0: 



(2.27) 



with he = h{rc) . 

These circular orbits, when they exist, have the minimum energy for a given 
angular momentum, and so are stable (at least within the present framework). If 
the energy is increased slightly, for fixed i, a broader class of orbits arises for which 
Pr ^ and so for which the radial coordinate oscillates between the nearest classical 
turning points, r±, defined as the two roots (on either side of Vc) of the condition 
£ = Vef[{r±) . For p > 5 the resulting orbits remain bounded away from r = since 
they oscillate around the radius of the stable circular orbit. For p < 5 hj contrast 
the bound orbits have radii which range from zero up to a turning point closer to 
r = than the position of the maximum of Ves ■ 

In general these orbits are not closed, since the period of the radial oscillations 
need not be commensurate with the angular period. In this case the orbits will be 
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ergodic in the sense that the system eventually passes through all points which are 
allowed by the given conserved energy and angular momentum. Since for p < 5 these 
orbits include radii near r = 0, the probe brane can and will pass near the source 
branes, and so is likely to be quickly annihilated. 

The existence of such unclosed orbits is not a surprise, given that similar orbits 
also arise for solutions to the relativistic Kepler problem (such as for the precession of 
the orbit of Mercury) . What we shall find which is a surprise is that for some brane 
systems the orbits are closed non-precessing ellipses, even once relativistic effects are 
included. 

2.3.2 Explicit Solutions 

We now turn to a more explicit characterization of the orbits. Notice that the above 
equations of motion imply that the derivative dr/dO takes the simple form: 

dr f pr 



(2.28) 



To solve for the trajectory followed by the particle we start from eq. ( p.28|) and 
use eq. ( p.22|) to eliminate pr in favour of the conserved e . Using the standard change 
of variable, m = 1/r, we have u' = du/dO = —pr/C- and so 



i2\ 



1/2 q 
~ h 



where now h(u) = 1 + ku'^. The orbit is then obtained by quadrature: 



uo 



dx 
F{x) 



where 



Fix) 



e' h'^{x) + 2eqh'^-'{x) + q' h 



s2 L<^-2, 



X 



- h'^-^nx) - fx^ 



1/2 



(2.29) 



(2.30) 



(2.31) 



2.3.3 Dp-branes 

Again 10 dimensional Dp-branes provide a particularly interesting special case, cor- 
responding to the choices rj = uj = 1 and q = q = ±1 • In this case the Lagrangian 
( p.l6|) becomes 

(2.32) 



L 



m 
J 



1 — /i (f2 + — q 



and so the canonical momenta become 



and 



pr = r 



1 + f h/r^^ 



1 — hr'^ 



(2.33) 
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The conserved energy then is 



l + h \ ^ + p. 



1)2 

_^ ,2 
2 ~^ " 



1/2 



q \ . (2.34) 



In this case eqs. (|2.3CI| ) and ( p.31|) defining the orbits become 



6-00= , , (2.35) 

Juo \/ A + Bx'^ -x'^ 

where A = {e"^ + 2 q e) / i"^ and B = ke^/C.'^ . What is striking about this solution 
is that it is precisely the same orbit as would be given by classical nonrelativistic 
motion in the presence of a central potential of the form Vc = —rrickc/r'^: 

I (r' + r'9') + h. (2.36) 



rric 2 V / 



The orbits for this Lagrangian ||12| are given by eq. ( |2.35| ), but with A = Ac = 2ec/tl 
and B = Be = 2kc/il, where as before Ec = Ec/rric and ic = Pe/rric are respectively 
the energy and angular momentum per unit mass. These correspond to replacing 
kc = 2k and taking the nonrelativistic limit: e = 2 + Ec with \ec\ <^ 1 . 

This shows that the orbits of the fully relativistic Dp-brane system have exactly 
the same shape as those of the nonrelativistic system — given by eq. (|2.36|) — with 



the corresponding orbital parameters read off by making the substitution {A^, B^) ^ 
{A,B). In particular, for bound states we take Sc < and so have ^4^. < 0. It 
follows that the bound Dp-brane orbits must also have A < 0. Now for Dp-branes 
the energy expression, eq. ( 2.34|) , implies e > for both q = ±1. Consequently 



A = [e"^ + 2 qe)/i'^ , can only be negative for antibranes, since then g = — 1. This is 
as expected since only if the probe is an antibrane does it experience an attractive 
net force. 

It is quite remarkable that the orbits for this system are exactly solvable at 
the same level as are those of the analogous nonrelativistic central force systems, 
particularly given that the Dp-brane Lagrangian includes an infinite number of terms 
when written in powers of 1/r and f^. We note in passing that although the shape 
of the orbits are unchanged from the nonrelativistic central-force problem, the time 
evolution of the branes along the orbits does change, as will be seen in more detail 
in a later section. 

2.3.4 Branes at Large-Distances 

To this point we have studied the orbits using the full solutions to the nonlinear 
field equations, which for Dp-branes can be done within a controlled approximation 
provided QsN ^ 1 . We now turn to the case where r is sufficiently large that only 
the weak-field form for the fields is required. This might be expected to apply in the 
case where N is not large. 
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For large r the virial theorem states that ~ k/r'^ and so if we expand in 
powers of k/r'^ we must also expand to the same order in powers of f^. Expanding 
the Lagrangian, eq. (|2.16|) , in this way gives (up to an additive constant) 




so we see that the leading dynamics is that of a particle moving in the presence of 
the central potential V = [q — l)k/r'^, which is attractive for q = —1. This is what 
is expected from explicit calculations using the string cylinder graph. The orbits for 
this case are given by eq. ( |2.35| ), but with A = and B = B^. 

We find in this way that eq. ( p.35| ) describes Dp-brane orbits for all r ^ . For 
QsN ^ 1 the parameters A and B are as given below eq. ( p.35| ), and for QgN < 1 
this becomes A ^ Ac and B B^. Although one might have expected the inclusion 
of the higher-order corrections in powers of v'^ and fc/r'' to more dramatically change 
the form of the orbits, this does not happen. Because the conditions rj = uo = = 1 
which ensure this simplicity for the brane/antibrane orbits (g = —1) also ensures the 
complete absence of an interaction potential between branes (g = +1), we conjecture 
that the independence of the orbit shapes from receiving relativistic corrections arises 
only for BPS branes and antibranes. A deeper understanding of when this occurs, 
and whether there is a connection with the BPS nature of the Dp-branes involved, 
would clearly be worthwhile. 

2.3.5 Closed orbits 

We next demonstrate the existence of brane orbits which close on themselves. Our 
interest here is not the circular closed orbits, which close on themselves in a trivial 
way and which arise for many kinds of branes. Rather we look here for the less trivial 
situation where the orbits vary in both the angular and radial directions, but these 
motions have commensurate periods. 

The existence of such orbits is suggested by the following two observations. First, 
we have seen that some relativistic branes (like the BPS Dp-branes in 10 dimensions) 
have precisely the same orbits as has the nonrelativistic central-force problem. Sec- 
ond, we know that motion in the nonrelativistic Coulomb potential, V oc 1/r, is 
along conic sections, with the bound orbits being ellipses. We focus therefore on the 
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case of source D6-branes and a probe D6-antibrane, for which d = 7 — p = 1 and 
q=-l. 



The probe-brane orbits are described by eq. (|2.35| ), which for = 1 may be 
integrated in closed form in terms of elementary functions. When A < the result 
is an ellipse 



a 1 - e" 
r = — — 

1 + e cos 6 ' 

where the orbit's semimajor axis, a, and eccentricity, e, are given by 

B ke 



(2.39) 



a = — 



2A 2 (2 - e) 



1 + 



4 A' 

S2 



1/2 



(2 - e) 



1/2 



(2.40) 



For comparison, using the nonrelativistic values A^ = '2Sc/(^1 ^^"^ = '^k^jl'l instead 
gives the usual expressions Oc = —kc/{2ec) and Cc = [1 + 2ecil/k'^]^^'^ . 

The time evolution along this orbit is obtained by integrating equation (2.14), 
giving 

{e + q)x + ek 



t-tr 



which can be integrated to give 



dx 



Qt = ip — e sin ip 



(2.41) 



(2.42) 



with the angle ip denoting the eccentric anomaly, defined by r = a(l — ecosip) or 
equivalently, 

a — r ke^ — 2 e{2 — e)r 

cosip = 



ae 



^k^e^-APe{2-e) ' 



The angles if) and 6 are explicitly related by tan | 
constant e which appears in (|2.42|) is: 



, ■ tan % . 

1— e 2 



e-1 
3-e 

The period of the orbit is given by: 

'3 -£ 



6-1 1^ U\2-e) 
3 — e " 



nn 27r 

T = — 



■nk 



2-e 



2-e 



TT J— (2a + 3A;) 



(2.43) 
Finally, the 

(2.44) 



(2.45) 



It is instructive to compare this with the nonrelativistic Kepler problem. In this 
case the time dependence is given by a relation similar to eq. ( ^^.421 ): 



fic^ = ~ sin (?/') 



(2.46) 
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where is related to 9 as before. The orbital period which results is Tc = 2tt /Qc = 



We see that the D6-brane and the nonrelativistic problems share the same func- 
tional form for both tlip) and r{6), but the parameters which appear in the two 
solutions differ. (For instance e ^ e for the D6-brane, while Cc = Cc nonrela- 
tivistically. ) This difference in parameters for the two solutions implies a different 
relationship is obtained between orbital features like a and e and the underlying 
physical quantities, like e and i and m. 

Although the orbits have the same shape, the differences in their parametric 
dependence on e and i also implies differences for how the period relates to the 
orbital shape and size. For instance neither Kepler's second nor third laws apply for 
D6-branes, since T^/a^ is not a universal quantity for all probe brane orbits, and 
since conservation of i does not imply that the product 9 is a constant along any 
one orbit. 

2.3.6 The Runge-Lenz vector 

Here we explore further the remarkable property that the closed orbits for the non- 
trivial D6-brane system, described by the lagrangian, ( p^.32D , are identical to the 
orbits for the standard nonrelativistic 1/r central force problem. 

For the nonrelativistic Kepler problem, the fact that the orbits are closed ellipses 
which do not precess is reflected by the existence of an extra conserved quantity: the 
Laplace-Runge-Lenz vector. Normally the existence of this vector is considered to be 
a special consequence of the 1/r potential, which is not shared by other central- force 
potentials. 

Nonetheless, despite the fact that the D6-brane Lagrangian contains many rel- 
ativistic corrections as well as a potential which is not purely proportional to 1/r, 
its orbits are also non-precessing ellipses. One might therefore wonder whether in 
this case there is also a generalised Runge-Lenz vector whose conservation ensures 
the absence of precession. We show in this section, by explicit construction, that the 
answer to this question is 'yes'. 

Since there are 3 spatial dimensions transverse to the D6 brane in 10 dimensions, 
it is convenient to phrase the analysis using the angular momentum vector, which is 
dual to the conserved antisymmetric tensor. Let us therefore consider the vector: 



where L = x A p is the conserved angular momentum vector and for the moment 
f{e) is an unspecified function of energy to be fixed below. The linear momentum 
is defined as the canonical quantity 




(2.47) 



P 



1 dL 
m 9x 



— hv 



1 



,2 ' 



(2.48) 
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and 



dL 1 



9x /iVl -hv'^ 



h V 



2 



1 gVl — hv"^ 



dh X 

Or r 



= £VrZ5Z|^- (2.49) 

2m or r 

where in the last equahty we used the expression for the energy eq. ( f^.22| ) as well as 

= 1. Using this expression to evaluate the quantity p in eq. (|2.47|) we obtain: 

1 dh 

A = 2fie)e'-^ [xA(xAx)]. (2.50) 
Now if we choose f{e) = and use the identity 

xA(xAx) = rrx-r'^x = "'"^ ^ ("") ' (2-51) 

we find 

A = -^r^^P). (2.52) 

2dr dt\r) ^ ' 

Now comes the main point. Only for the particular case h = 1 + k/r is the 

right-hand-side of this last equation a total derivative, and so can be used to make 

a conserved quantity — the generalised Runge-Lenz vector: 

1 A; X , , 

A = -pAL---. 2.53 

e"^ 2 r 



A-x = Arcos6' = — - — (2.54) 



Notice that, as expected, A ■ L = 0, and 

£^ ~2 

which exactly reproduces the elliptic trajectory r{6) as before, and sets A = ke/2. 

As in the nonrelativistic case, the Runge-Lenz vector indicates a further symme- 
try beyond rotation invariance, and because it lives within the plane of the ellipse, 
pointing along the semimajor axis, its conservation forbids the orbits from precessing. 

2.4 Quantum Dynamics 

We next describe the quantum behaviour which follows from the lagrangian, ( |2.32| ). 
We do so for its intrinsic interest, even though all of our later brane applications will 
be well within the classical domain. 

Recall that for Dp-branes the Hamiltonian is the following function of the canon- 
ical momentum ^ 

Ti. = — {^\/rri?^Ph^ — mq^ . (2.55) 

The Schrodinger equation for this Hamiltonian, Ti^ = E^^ , therefore becomes 

^ i^-^m? + hp'^ -mq^ ^! = . (2.56) 
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This may be written in a more useful form by isolating and squaring the square root 
on the left-hand-side, leading to 

{m^ + hp^) ^ = {Eh + mqf ^. (2.57) 

Because we have squared Ti in what follows we will find solutions having both signs 
for E. Since Ti is manifestly positive for g = ±1, we must drop all solutions for 
which < (notice that the time-dependent Schrodinger equation will include two 
time derivatives indicating positive and negative frequencies). 

For 0^ = 1 and h = 1 + k/r'^ this can be further manipulated to be written as: 

^2 ^\ 

J- j ^ = E{2mq + E)^ . (2.58) 

The utility of this last form comes from its similarity with the corresponding nonrel- 
ativistic central-force problem: 

* = (2.59) 

which becomes identical to eq. (|2.58|) after the replacements /i —>■ 1/2, g"^ E'^k, 
and Ec^ E{E + 2qm) . 

In particular, for D6 branes we have d = 1, and so eq. ( |2.58| ) has the same form 
as the Schrodinger equation for the Hydrogen atom, for which the solutions have been 
known since antiquity. In this case the energy eigenvalues consist of a continuum, 
Ec > 0, plus a discrete spectrum of bound states labelled by the quantum number 
n=l,2,...: E, = -g^^I/{2n')^ 

Making the substitutions which convert this spectrum to the D6-brane problem, 
we find the following energy spectrum: 

E{E + 2qm) > , (continuum) (2.60) 

E'^k'^ 

E{E + 2qm) = — ' (discrete). 

Requiring an energy spectrum which is bounded from below implies the continuum 
branch becomes E > if q = 1 and E > 2m if q = —1. 

The energies for the discrete spectrum can be explicitly solved to give the fol- 
lowing expression in terms of the level n and the constant k: 



Ein) 



A;2 



1/3 / I \ 1/3" 



' m? -n? qm \ I m? 'n? qm 



(2.61) 

Notice that the square bracket in this expression is negative if g = +1 and is positive 
if g = — 1 . Consequently for g = +1 these discrete eigenvalues are spurious and must 
be discarded. For g = — 1 the discrete eigenvalues are positive, and approach the 
continuum {E = 2m) a.s n —>■ oo. The first 20 energy levels which follow from 
eq. (|2.61|) are plotted against n in Figure |^. 

^Recall our use of fundamental units, for which h = c = 1. 
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The characteristic size of the 
corresponding discrete energy eigen- 
states in the transverse dimensions 
is set by the D6-brane analog of the 
Bohr radius: as = l/{g^ After 
the appropriate substitutions this 
becomes 



as 



(2.62) 



Energy 



For subsequent apphcations it 
is interesting to see how these ex- 
pressions behave in the decompact- 
ification hmit, for which the spa- 
tial world- volume, Vq, of the brane 
becomes very large. Because m = 
TqVq, this corresponds to the large- 
m limit, with n and the interaction 
constant k = c^QsNlg held fixed. 
We see that in this limit the energy 
and the Bohr radius of the bound 
states scales like 



0.66 



0.65 



0.64 



0.63- 



0.62- 



0.61 



0.6 



8 10 12 14 16 ie 



20 



Figure 2: The energy of the first 20 bound states 
of the branonium. We have taken m = 1/3 and 



E oc 



m 



1/3 



^1/3 



p/3 ' • ^2.63) 

We see from this that the quantum bound states are not reliable in the limit of large 
Vq , for fixed N , because in this limit they imply as 

2.5 Compactification 

We have found potentially long-lived orbits for probe antibranes moving in the field 
of source branes, provided d = 1 and so they move in 3 transverse dimensions. 
For 10 dimensions overall, this is what singles out D6-branes for special attention. 
For later applications to brane- world cosmology we will be interested in branes whose 
world-volume is 4 dimensional, and so we pause here to consider how this may be 
achieved from the above solutions by compactification. 

2.5.1 Toroidal Reduction 

Our goal in this section is to dimensionally reduce the action, eq. (|2.1| ), from D 
to D' = D — 1 dimensions. The point is to present a consistent lower-dimensional 
truncation of this action, in the sense that any solution of the lower-dimensional 
equations of motion are also solutions of the equations of the full system, eq. ( |2.1| ), 
following the construction of ref. [ITU . 
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To this end imagine that one of the coordinate directions, denoted z , is periodic. 
In this case we write the D- dimensional metric, in terms of a D'-dimensional 
metric, ds\,t plus a scalar field, y?, and a vector field, Bm as follows: 



ds\^ = e'^^ dsl, + e^^'P {dz + Budx^^f , (2.64) 

with h = —{D' — 2) d and = 1/[2{D' — 1)(-D' — 2)] chosen to ensure canonical 
Einstein and ip kinetic terms. The D-dimensional {n — l)-form field, A[n-i], is 
similarly decomposed into D'-dimensional {n — 1)- and {n — 2) -form fields, -B[n-i] 
and -B[„_2] , according to 

A[n-i] = B[n-i] + B[n-2] A dz . (2.65) 
With these choices the action, eq. ( p.l|) , dimensionally reduces to 

Sd' = - j d^'x [l Ro^ + \ {dct>f + i {d^f + \ e-2(^'-i)»^ 



I ^ -2(n-l)a(/3+a(^^'2 2{D'-n)a95+a ^2 



(2.66) 



where the field strengths appearing in this expression are defined by 

B[2] = dB, G[n-l] = dB[n-2], G'^n] = dB[n-l] — G[n-l] A B . (2.67) 

This lower-dimensional action may be consistently truncated in three distinct 
ways, in the above sense of consistency. The three ways are obtained by setting to 
zero all fields except the D'-dimensional metric and one of the following pairs of 
form-fields and canonically- normalized scalars ||1CI|| : 

Option A: B and (f)[ = ip 

^/4:{D' - n)2 + 

n n D A A' -2in - 1) d + a(j) 

Option O: -D[n-i] and (p^-i = 



/4(n-l)2a2 + a2 

Each of these reductions produces a -dimensional action which again has the 
form of eq. (|2.1|) , but with parameter a = a^' now given in terms of the higher- 
dimensional parameter by: 

Option A: = —2{D' — l)aD implying A' = 4, 



Option B: a^' = y 4(/^' - nyd"^ + al implying A' = A , (2.69) 
Option C: az?' = ^ 4(n — lyd"^ + aj^ implying A' = A , . 

The two main choices are now whether to apply this reduction to the case where 
the compactified dimensions are transverse or parallel to the source and probe branes' 
world volumes, and so we now consider both of these options. 
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2.5.2 Parallel Dimensions 

The simplest option is to compactify dimensions which are parallel to the branes' 
world volumes by wrapping the branes about a torus. Wrapping D6-branes about a 
3-torus of radius R in this way would produce a brane world which is effectively 4 
dimensional when viewed on scales large compared with R by virtue of the standard 
Kaluza-Klein mechanism. Such a toroidal compactification is straightforward since 
the solutions we have considered are invariant under translations along the directions 
parallel to the branes. In particular, it is clear that the solution, eq. (|2.3| ), remains a 
solution to the full higher- dimensional field equations if some of the coordinates 
are required to be periodic. 

Writing = z and directly comparing eqs. ( p.3|) with the Kaluza-Klein ansatze, 
eqs. ( 2.64| ) and (|2.65|) , allows the field configurations to be written purely in terms 



of lower- dimensional fields as follows: 

dsl, = h-^' dxl, + h^' dy\ e^ = h\ e'^ = /^^ (S[„^2])oi...p-i = C(l - /^"') , 

(2.70) 

where k and C, are as before (c./. eq. (|2.4| )) and 



^ 7' = 7-2ap, 7' = 7 + 2ap. (2.71) 



2{D' -2)a 



Notice that although this solution resembles truncation option B of eq. (|2.68|) , it 
is not simply a special case of the solutions we have considered so far, with the 
replacements D ^ D' , p p' and a — a'. It is not because in this solution 
both of the scalar fields and vary nontrivially in the transverse dimensions. (In 
particular, the combination orthogonal to 0^_2 of eq. ( p.68|) does not vanish.) 



By starting from a given p-brane solution in D dimensions and successively com- 
pactifying k of the parallel dimensions we obtain new configurations which describe 
p — k branes moving in D — k dimensions, without changing the number of trans- 
verse dimensions and so also not changing d. (These are the 'diagonal' reductions 
of ref. ||lOl-) By construction this reduction leaves the values of the parameters r], u 



or C which appear in the probe brane action unchanged. 

Reducing the BPS Dp-branes in 10 dimensions in this way then leads to a family 
of p — k branes moving in 10 — A; dimensions, for which rj = u = ( = 1 , and so for 
which the orbital properties are particularly simple. The choices p = 6 and k = 3 
give the particularly interesting case of 3-branes moving in 3 transverse dimensions. 
For these d = 1 and so the orbits are the closed ellipses described earlier. 

2.5.3 Transverse Dimensions 

To this point we have treated the transverse dimensions as if they were infinite 
in extent. Although this type of geometry might have brane- world applications if 
gravitons can be confined to our brane (perhaps through the warping of spacetime). 
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we need not rely on this when exploring the various implications of branonium. If 
not, we must also address the dimensional reduction of the transverse dimensions. 

Reduction of the transverse dimensions is trickier because the field configurations 
in these directions are not translationally invariant. For BPS branes, such as the 10- 
dimensional Dp-branes, this can be dealt with using a generalization of the solutions 
we have considered to source branes which are not all located at the same point. This 
more general solution is obtained using the same expressions, eqs. ( p.3|) , as before 
but with the harmonic function of eq. (p.6|) replaced by the new expression 



h{y) = l + Y. ^ , (2.72) 

where |y — y^P = Smniv — yi)"^{y — UiY ■ Given this solution a toroidal compact- 
ification of the transverse dimensions may be performed simply by including the 
appropriate image charges for the source branes, which serve to ensure that the 
transverse dimensions become periodic [|1^ . 



There are two cases of this kind of reduction which concern us, depending on the 
relative size of the compact dimensions and the size of the brane orbits which are of 
interest. If some of the transverse dimensions are very small compared to the orbital 
sizes of interest, then we may integrate out these directions and work in the limit 
where their radius is taken to zero. If k transverse directions are so reduced it has 
the effect of reducing D ^ D' = D — k without changing d , and so gives new brane 
solutions for which d d' = d — k . (These are the 'vertical' reductions of ref. ||10|| .) 

For the brane orbits of previous sections we must take d > 2, however, and 
so at least some of the transverse dimensions must have radii, R, which are larger 
than the orbital sizes, a, of interest. In this case the orbital analysis of the previous 
sections must be considered to be only approximately valid, up to corrections which 
are of order {a/RY for a <^ R. For a ^ R the influence of the image sources can 
significantly alter the probe-brane motion. In our later applications we imagine at 
least three of the transverse dimensions to be stabilized with sufficiently large radii. 



3. Stability 

We now turn to the stability of the orbits described in the previous section (see [IT3 



for a related discussion). There are several kinds of potential instability to consider. 
To this point we have treated the probe-brane motion as if the branes were perfectly 
rigid, with only their centres of mass free to move. The first instability issue to 
address is the extent to which this is a good approximation, inasmuch as other kinds 
of brane motion might become excited by the movement of the centre of mass. To 
address this at the classical level we look for tachyonic modes among the linearized 
fiuctuations of the brane position, y{^,t), about the field configurations of interest. 
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We argue that this kind of instabihty need not be present so long as the source and 
probe branes are kept at distances which are much larger then the string length. 

The second stability issue asks 
whether the brane motion produces 
radiation of massless bulk and/or 
brane fields, thereby radiating away 
energy and so causing the brane or- 
bits to decay. This kind of instabil- 
ity does occur, and we here estimate 
how long it takes to destroy a given 
brane orbit. 

All told we identify three sepa- Figure 3: Nucleation of a throat between a 
rate kinds of instabihties to which brane and and antibrane. 
orbiting branes are subject. 




• Brane annihilation, due to the perturbative open string tachyon. If the brane 
and antibrane get within a string length or so of one another, then they can an- 
nihilate with the antibrane being absorbed by the source branes, ultimately 
leaving a configuration of — 1 branes. The energy released is then radiated 
into modes in the bulk and on the remaining source branes. We do not pursue 
this in detail because our present interest is in branes which are sufficiently 
well separated as to avoid this. Notice that, depending on their number and 
size, for compact transverse dimensions sufficient separation of the branes can 
require them not to be too differently wound around the extra dimensions. 



A non-perturbative brane decay is possible, a la Callan-Maldacena 0]. These 
are solutions to the Euclidean equations of motion. They can be interpreted 
as sphalerons associated to the decay of a brane- antibrane system. A throat is 
nucleated that joins the brane- antibrane pair. Then it expands and transforms 
the brane world-volume into kinetic energy of the expansion front (see figure 
^.5.31 ). We do not pursue this instability further since its nucleation probability 
is proportional to e'^l^" , and so is exponentially small for small string coupling. 



• Radiation into bulk modes, such as the graviton and dilaton, or brane modes 
can occur and arises at low orders in the string coupling. We expect this to be 
the dominant instability for large interbrane separations and weak couplings, 
and estimate the time for orbital decay below. 



We will then establish first the stability of the system against bending and con- 
centrate later only on the instability due to radiation. 




Figure 4: The fields y^{(,) on the brane parameterise its fluctuations. Instabilities corre- 
sponding to brane bending or the classical orbits can be analyzed from the mass of these 
modes in the quadratic approximation. 



3.1 Stability Against Brane Bending 

Instead of considering only explicitly straight and parallel branes, imagine allowing 
them to bend so that the probe-brane position relative to the centre-of mass is 
y\t,^^) = y^ — Vo, where as before {t,^^} are the coordinates along the brane (see 
figure Nontrivial dependence of on corresponds to brane bending, while 
perturbations depending only on time, i.e. := y\t) , describe the rigid motion of 
the brane centre of mass. 

A stability analysis is performed by expanding the action about a particular 
background solution, = yo(t) + 6y\^,t), with the background yo(t) taken to 
describe the rigid brane motion described in the previous sections. The expansion of 



L about this solution defines the system's Routhian [12|, which to quadratic order is 



-R = dM W + ^y' W + ^y'^y' 



where 



and 



CrAyoit)] 



dy ^ dy^ 

A'^-iyoit)] 



d^R 



d^R 



d{d^y') dyi 



d{df,y') d{dyy^) 



(3.1) 



(3.2) 



(3.3) 



One must analyse this effective Lagrangian looking for unstable perturbations, given 
the known functions for the quantities A'^^ , B^^ and Cij . 
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For the present instance, we may for simplicity focus on the large-separation 
limit and so take the non-relativistic approximation 



m 



(3.4) 



The stability analysis greatly simplifies in this case because simplicity of the deriva- 
tive terms implies the quantities A'f^J and B^- are independent of fields, and so are 
also independent of time. As a result it also follows that B^ - = 0. 

The simplicity of the kinetic terms allows a further simplification, since it implies 
gradients contribute a strictly positive amount to the energy density. This means 
that it is the zero-momentum modes in a Fourier decomposition, 6y^{k) oc e**^'^, which 
are the most unstable. It follows that if these zero-momentum modes are stable, then 
so must be all of the modes having nonzero momentum. This stability of nonzero 
momentum modes has a simple physical interpretation. On one hand, if the probe 
brane were to bend slightly towards the source branes, then tidal forces from the 
source branes will act to increase the amplitude of this bending. The amount of 
this tidal disturbance will fall with increasing brane separation. On the other hand, 
the brane tension resists all such bending and easily wins so long as the interbrane 
separation is large compared with the string length. 

We now specialize to the stability analysis for the zero-momentum modes in a 
purely time-dependent background. The effective Lagrangian for the perturbations 
dy^ = — Vq then simplifies to 



L 



-R = ^AjSyW + B.,5y%^ + l^.SfSy^ 



(3.5) 



where C,- 



grangian, eq. 



yo 



dy^ dy^ 



and A 



d^R 



yo 



V 



dij^dy^ 



Applied to the la- 



yo 



, we get the following Routhian expansion for the perturbations: 



m 



2 



(Srf + rl {Sef + {5y 



i\2 



iv - q) 



kd{d + l] 

1^+2 

'o 



3 I ^ 



{5rf. 
(3.6) 

Only the radial mode appears nontrivially in these manipulations, leading to the 
'mass': 

M-,=-m{^-q)^-^ + -. (3.7) 

Our interest is in — g > 0, such as occurs for probe Dp-antibranes, and so the first 
term in is negative while the second term is positive. Which term dominates 
at large r depends on the relative size of the power of r in the two terms, and so 
< for large r when d < 2 and > for large r when d > 2. In the 
intermediate case d = 2 the sign of at large r is the same as the sign of the 
combination — mk{ri — q) d{d + 1) . For the Dp-branes this means that we have 
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classical instability at large r for p > 5 and classical stability for p = Q. This is as 
expected since we know that deviations from the classical orbit in the Kepler problem 
for (i = 1 are purely oscillatory. 

Since D6-branes can be supersymmetric, the oscillations about the stable con- 
figuration fill out an entire supermultiplet of the appropriate unbroken supersymme- 
tries. Seen from 4 dimensions on the brane, this could be an = 4 SYM multiplet, 
depending on how many supersjTiimetries are broken by the compactification. 

A similar analysis can be done to study the stability of the brane/antibrane 
system against developing a relative rotation as their centres of mass move. This 
does not occur because parallel branes furnish a local minimum of the potential 
energy as a function of angle fill]. For branes embedded in infinite lOD spacetime, 
the oscillation time scale for small rotational perturbations can be quite large, due 
to the necessity to rotate an extremely long brane. In this case one should check 
how this oscillation time scale compares with the orbital period calculated above. 
For real applications this time scale is less important, however, because the branes 
then move within compact transverse dimensions. For toroidal compactifications, 
for instance, branes cannot rigidly rotate relative to one another without a large 
energy cost. This cost arises because local rotation of the branes necessarily requires 
the branes to bend somewhere, due to the necessity that they cannot change their 
windings about the cycles of the torus. 

We do not expect this result to change when relativistic corrections are taking 
into account, and so we expect the relativistic system to also be stable (for d = 1) 
provided the interbrane separation is kept well away from the string length. We now 
turn to the orbital energy loss due to radiation of various string modes. 

3.2 Radiation into the Bulk 

The most important decay mode for large orbits at weak coupling is energy loss 
through radiation into the various massless fields in the problem. For the classical 
motion which we are considering, this radiation causes the orbits to slowly shrink in 
size, much as the binary-pulsar orbit has been observed to do due to gravitational 
radiation. Once the orbital size reaches the string scale the much quicker open-string- 
tachyon instability takes over, causing a catastrophic annihilation of the antibrane 
with the brane. Our purpose in this section is to estimate the time which is required 
before this catastrophe takes place. 

3.2.1 Radiation Rates 

Unlike the situation for the binary pulsar, the moving probe antibrane has several 
types of charges which couple its motion to fields whose radiation could carry away 
energy, including at the very least the bulk gravitational, {n — l)-form and dilatonic 
fields. For nonrelativistic motion, when all other things are equal, the main carriers 
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R + 



of radiated energy are the lower-spin fields [17|.^ This is because in this limit spin-0, 
-1 and -2 particles respectively dominantly couple to the orbiting brane's monopole-, 
dipole- and quadrupole-moments, with each higher multipole suppressed by addi- 
tional small factors of powers of 1/r. We therefore base our estimate on radiation 
into a bulk scalar field like 0. 

Suppose there are D dimensions which are large compared with the brane orbits, 
of which d are parallel to the various branes. The case of special interest is d = 4 
and d = 1, which corresponds to D = 7, but for the present purposes we need not 
restrict ourselves to this case. We imagine the probe brane orbit to be circular, with 
orbital radius r and orbital frequency Q. 

Recall that (p appears in the brane action through the overall factor e~'^m, 
with (as usual) m = TpVp, where the spatial volume of the brane enters because of 
the overall translation invariance within the brane dimensions. This action is given 
in units where the scalar and Einstein kinetic terms are S = — ^ J d^x ^J—g 

((90)^ , and so we see that the 0-brane coupling is proportional to K^-m, where 
= SvrG^j/V^ relates to the D -dimensional Newton's constant. Furthermore, 
since for slowly-moving objects Lorentz-invariance precludes any radiation unless the 
source accelerates, the 0-emission amplitude should be proportional to the brane 
acceleration, a^- 

Keeping track of dimensions and factors of 27r, the power radiated into scalars 
per unit brane volume therefore has the schematic form 

P^27T4m'alJ 6{E-n)j^^^^^^, (3.8) 

where the final integral is over the phase space of the emitted scalar. The Dirac 
(5 -function is meant to express the fact that energy conservation demands the domi- 
nant energy of the outgoing scalars to be set by the dominant frequency, Q, of the 
time-dependent source. It is this energy which sets the scale of the phase space 
integrations. The power of E in the integral is chosen on dimensional grounds. 

We estimate the phase-space integration by replacing it with the appropriate 
power of the relevant energy scale, Q, and by counting the factors of 27r which are 
produced by the angular integrations in the transverse dimensions. This leads to the 
estimate that the phase-space integral is of order (27r) t*^'^"'^^/^] , where [k/2] 

denotes the integer part of k/2. We find in this way the following estimate for the 
total power radiated: 

Pr^4m'al (%^] , (3.9) 



(2vr) 



with V = D -d-1- [i(D - d)] . 



^This observation also applies to the binary pulsar in the presence of very light scalar and vector 
particles, as may be seen, for example, in ref. [na. 
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For three transverse dimensions we would have D — d = 3 and so [j{D — d)] = 1. 
This would be appropriate for point particles (0-branes) moving in 4 dimensions, in 
which case we would have D = 4 and d = 1 and kt = 1^4 is the ordinary 4D 
gravitational coupling. The above estimate then reproduces the results of more 
careful calculations, which give: 

P = ^ (K4m)\l . (3.10) 

OTT 

This same formula applies equally well to the case of more general branes orbiting 
in 3 transverse dimensions, provided we use the appropriate choice for kt 

P^l-^m'al =j-nlT^alV,. (3.11) 

Notice that in all cases P scales with the spatial volume of the orbiting branes, 
as is required by their translational symmetry, and so the power radiated per unit 
energy is independent of V^. As a consequence of this Vp drops out of the time 
required for an orbit to decay, as we shall shortly sec in more detail. Also of note 
is the dramatic phase-space suppression (by powers of Qls) which arises when there 
are more than 3 transverse dimensions. 

3.2.2 Orbital Decay Time 

We now turn to the time required for a circular orbit of radius r to decay, given the 
above estimate for the power radiated. We focus on the case d — 1, since it is for 
this case that these stable orbits are most likely to exist. For the present purposes we 
write the power radiated as P = K al , with K containing the various dimension- 
dependent constants. 

To evaluate the acceleration it suffices to work in the limit of large orbits and 
small velocities, where the brane Lagrangian becomes (up to an additive constant) 

L^^-Vir) (3.12) 
where V — —2km/r. The acceleration is then 



]_(dV_V 4/^2 
\ dr 



and so the power radiated becomes 

P^-—^ Km^ — = — -— . 3.14) 

This last equality uses the nonrelativistic expression which relates the energy and 
radius of a circular orbit: E — —2mk/{2r) — —mk/r. 
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The time taken for the radius to change from to is therefore 



"^f dE m'^k^ r^^dE 



3 ^3^ 



r,; — r 



Taking r f Ig ^ ri and comparing r with the initial orbital period, T = 27r a/ rf /2k 
gives an estimate of the number of orbits traversed before decay: 

3/2 

A/" = - = ^ . (3.16) 

T 12tx V2kmK 



To see how these results scale with the microscopic quantities, we may use the 
D6-brane results m = TqVq = Vq/ {{27ry I gs) and k = c^QsNlg = gsNls/2 and 
K = kIq/(8itVq) with kIq = (2n)^TTg'^ls^ . Using these results the number of orbits 
before decay becomes 

2 fvA'/' 1 



Ar=— / — . (3.17) 



This dependence is relatively easy to understand. 



• As already discussed, Vq does not appear in eq. (|3.17|) because P and E are 
both proportional to the brane volume. 

• The string coupling dependence arises from several sources. First notice that 
the orbital energy does not depend on g^ because the contribution cancels be- 
tween the interaction constant k and the tension Tp. For a similar reason gg 
drops out of the combination m^/t^ in the radiation rate. The entire contri- 
bution of gs to the radiation rate therefore arises through the acceleration, 
which vanishes if — 0. Taken all together this gives a l/g"^ dependence to 
the lifetime, which combines with the g^^ arising from the k^^"^ in the orbital 
period, to make N (x gs ^''^ . 

• The factor is the relative charge of the set of branes to the antibrane. The 
bigger is iV, the stronger is the acceleration of the probe brane and so the 
shorter is the branonium lifetime. 

We see that the branonium lifetime is maximized the smaller we take g^^ N and 
the larger we take the orbit to be. As discussed in previous sections, if we wish to use 
the full nonlinear field equations rather than just their linearized form, consistency 
requires ggN ^ 1, and so a long-lived orbit would require gs and r/ls to be chosen 
that much larger. 

Of course if the brane were in the quantum limit, this orbital decay instability 
could be prevented just as is done for electromagnetic decay in atoms. Unfortunately 
this mechanism is not available to us for branes in the limit of large Vp , since for these 
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the stabilization occurs at distances of order the Bohr radius from the source branes. 
As we have seen earher, in the decompactification hmit the stabihzing bound states 
he beyond the domain of vahdity of our approximations, since their size is smaller 
than the string length. It might nonetheless be of interest in some contexts, such as 
in the very early universe if we assume that all dimensions start with a size close to 
the string scale. 

3.3 Radiation into Brane Modes 

We have seen that probe brane orbits decay due to radiation into bulk modes, and 
but this radiation could also be accompanied by radiation onto the brane. We next 
perform an estimate of this rate, along the lines of the one done earlier for bulk 
radiation, with the conclusion that radiation into brane modes is subdominant to 
bulk radiation. 

We use for these purposes the coupling of the brane-bound gauge field, given by 
the (String-Frame) Born-Infeld action: 

Sbi = -TpJ d^i e'^ ^ - det(W) [l + + •••], (3.18) 

where ^f^^ = QMNd^x'^d^x^ is the brane's induced metric and J^^j, is the gauge field 
strength. We see that the relevant coupling is in this case kb = Tpls^e-'t'\ where 
(pb = (f){yh) is the dilaton field evaluated at the position of the brane. 

Although radiation into the bulk requires non-vanishing acceleration, this is not 
sufficient to produce radiation into brane modes. This is because brane radiation 
arises due to the time dependence of the fields which probe-brane-bound observers 
see as a result of their motion through the fields of source branes (much as occurs in 
mirage cosmology). Since this time-dependence only arises if the interbrane distance, 
r, changes with time, brane radiation does not occur for circular brane orbits. 

The role of al in the bulk radiation rate is in this case therefore played by 
(dh/dt)^ ~ {2kr/r'^Y ~ {2kVte/rY ^ where e is the orbital eccentricity (which we 
take to be small). Keeping in mind that translation invariance requires the power 
radiated to be proportional to Vp and that the particles be produced in pairs — with 
opposite momenta to conserve momentum — we obtain the following estimate for 
the power radiated per unit spatial volume 

where the power of energy is included on dimensional grounds and V = 2p — 1 — [p/2]. 

Choosing 3-branes, and using kb = T^h^ — ^/gb, where gb = gs^'^'' is the 
string coupling at the position of the brane, gives \p/2] — 1 and so 

(3.20) 
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Here we use the large-r results i7 = ^J2k/r^ , gb = gs[l — {2kK/r) + ■■■], as well 
as the relation 2 ~ ggNls- The dependence of this estimate on Q, e, and k/r is 
borne out by a toy calculation of the rate of change of the vacuum energy density on 
the brane which is given in the appendix. 

Several conclusions may be drawn from this estimate: 



Comparing the power radiated into the bulk ( |3.11| ) and ( p.2q ) gives the ratio 



of power radiated into the bulk and to the brane. For D3-branes moving in 3 
transverse dimensions this ratio becomes: 



Pbuik ^7K2vr)3«:? T| al \dt J gli2T:f \ r J {2TiYgl \ r , 

(3.21) 

This shows that the greater part of the emitted radiation goes into the bulk. 

Although most of the radiation goes into the bulk, particles are produced on 
the branes, and this has possible implications for cosmology. Provided the 
number of source branes differs from the number of probe branes, some branes 
survive even after the orbit eventually decays. If, for instance, we imagine 
ourselves being trapped on one of these remaining branes, these particles can 
equilibrate with observable particles at later times. Depending on the details, 
this could provide a source of reheating and/or baryogenesis. For a branonium- 
inflation scenario these produced particles are likely to be inflated away as they 
are produced, and so are unlikely to obviate the necessity for other sources of 
reheating, such as due to the brane-antibrane annihilation itself. 



4. Conclusions 

We have seen that the bound states of brane systems have very interesting properties. 
The present work should be considered only as the beginning of their exploration. 
Within a brane world context we find it appealing to imagine our entire universe as 
being, albeit briefly, a member of a 'solar system' of universes. 

Our analysis has included only branes and antibranes which share the same 
dimensionality and are parallel, which makes a closer analogy with positronium than 
with a stable atom. It is straightforward to generalise our discussion to a system of 
branes having different dimensionality and/or branes at angles, for which the effective 
interactions can still be attractive. We would expect the behaviour and instability to 
annihilation of these systems to be very similar, although a detailed analysis is left 
for the future. Similarly one might extend this idea to possible bound states between 
branes and orientifold planes. Other systems, including lifts to 11 dimensions and/or 
T- and S- dual versions might also be considered. 

Many open questions remain. One is a more fundamental understanding of why 
D-brane/antibrane orbits do not precess, having the same shape for fully relativistic 
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systems as they have in the non-relativistic hmit. There may exist a mapping between 
the two Lagrangians in a similar way as recently found in ref. for the case of a 
massive particle near the horizon of an extreme Reisner-Nordstrom black hole. Notice 
that Bertrand's theorem |12| states that the only systems giving rise to closed orbits 
in the central force problem are the 1/r potential and the harmonic oscillator. The 
D-brane systems we consider here are not counterexamples to this theorem because of 
the canonical kinetic terms which the proof of Bertrand's theorem assumes. It would 
be interesting to see how other classical results, such as the virial theorem, the three- 
(or more) body problem or chaotic orbits, generalize in the present instance. 

Our quantum treatment of the branonium system was very basic. A better un- 
derstanding of how to properly describe the system at the quantum level and the 
physical interpretation of the corresponding wave function would be desirable, per- 
haps with implications to quantum cosmology. It would also be interesting to know 
in more detail the back reaction of the probe brane on the bulk metric, both to better 
understand the = 1 brane/antibrane problem and to have an explicit microscopic 
derivation of the long-distance fields which are of interest for cosmological applica- 
tions. Because the branes behave as point masses in the transverse dimensions it 
is clear that for 3 or more transverse dimensions the bulk gravitational (and other) 
fields in the extra dimensions can be made negligible simply by keeping them well 
separated. Experience with co-dimension 1 branes can be misleading in this case, be- 
cause the bulk fields which these give rise to can instead be homogeneous spacetimes 
like de Sitter or anti-de Sitter space. 

Finally, the applications to cosmology are potentially many. Inflation can be 
analyzed in a way similar to ref. [Q. In this case the system has several novel 
features. For instance it provides an unusual slow roll because it is the angular 
variable, 9, which rolls even though the potential is exactly flat in this direction. 
Nonetheless, the constant inflationary energy density can be expected to evolve as 
the orbital radius changes due to energy loss into bulk radiation. This changes the 
analysis in an interesting way, since r is not related to dV/ dr in the usual expression. 
Furthermore, the fact that the orbit is unstable and decays is likely to make the brane 
motion appear as a very unusual 'fluid' from the cosmological point of view. As for 
other forms of brane/antibrane inflation, the end of inflation would be occur once the 
antibrane falls into the source branes and annihilates, providing the usual extremely 
rich realization of hybrid inflation. 

The orbital motion of the probe brane might also provide other novel cosmolo- 
gies, such as a particularly interesting example of a mirage cosmology. For elliptical 
orbits the periodic increase and decrease of the interbrane separation gives rise to 
a cyclic cosmology for the trapped brane observers. Because the small size of the 
extra dimensions implies the period of these cycles is microscopically short for stan- 
dard compactifications, this is likely to be relevant only for very early cosmology. 
The particle production which is produced by the time dependence of the induced 
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fields on the branes in these cosmologies may have implications, perhaps providing a 
geometric realization of the Affleck-Dine mechanism for baryogenesis. Because these 
speculations take us considerably beyond the scope of the present paper, we devote 



a companion article |jT9| to a more detailed discussion of these cosmological issues. 



5. Appendix: Quantum Energy on the Brane 

We here briefly record a toy calculation of the time dependence of the quantum stress 
energy density on the brane which arises due to the time dependence of the induced 
brane metric. 

Our toy model consists of a conformally coupled scalar, in 4 spacetime dimen- 
sions, defined by the action 

C 



--\{dvf~\R^\ (5.1) 



We evaluate the quantum stress energy tensor for this scalar using a time-dependent 
but conformally-flat metric 

ds^ = h-^'^{-de + di'') , (5.2) 

where h = 1+k/r and r = tq 1+e cos{Qt) . Here e and Q represent the eccentricity 
and the angular frequency of the radial oscillations for a slightly eccentric elliptical 
orbit. 

Now, for any theory which is conformally trivial, the expectation value of the 
stress energy tensor, T^j, , is determined by integrating the conformal anomaly, lead- 
ing to the expression 

b 



14407r2 



(^R■,^lv — ag^iiyOR + ( « — ^ ) RRfiu 



+ ( ~4 8 J 9lJ.uR^ + ^RfiaRu - -^QnuRapR"'^ 



(5.3) 



Here the constants a and h depend on the spin of the field in question. For a 
conformal scalar fields we have a = — | and & = 1 , while for a massless vector field 
we have, a = — 3 and 6 = 62. 

Inserting the metric, eq. (p.2|), into eq. (^.3|), making use of the relations uj = 



^y2kjr^ , and m = TpVp we find (to leading order in I/tq and e): 

P = (Too) = ^^^^ [2a + {b- 3a) cos'm] . (5.4) 

An estimate for the energy production rate on the brane is then 
P dp Q^k'^ 



V3 dt 12807r2r2 



[-b + 3a) cos(fit) sin(fit) . (5.5) 



Notice that this scales like {2k / r)^ , just as does our estimate in the main text 
for the energy deposition rate for the brane. 
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